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Motivated by questions of present interest in nuclear and condensed matter physics we consider 

^^ , the superposition of a diagonal matrix with independent random entries and a GUE. The relative 

0\ ' strength of the two contributions is determined by a parameter A suitably defined on the unfolded 

scale. Using results for the spectral two-point correlator of this model obtained in the framework of 

the supersymmetry method we focus attention on two different regimes. For A <C 1 the correlations 

Tn ' are given by Dawson's integral while for A S> 1 we derive a novel analytical formula for the two-point 

fT , \ function. In both cases the energy scales, in units of the mean level spacing, at which deviations 

from pure GUE behavior become noticable can be identified. We also derive an exact expansion of 

the local level density for finite level number. 
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I. INTRODUCTION 



Transitions from regular to chaotic fluctuation properties have been observed, experimentally and theoretically, 

in many different areas of physics. As examples, we mention topics in atomic (see, e.g. Ref. nl) and nuclear (see, 

e.g. Refs. g-§) physics, the relation between level statistics and localization in condensed matter physics, and, of 

'"^ , course, quantum chaology. Although the physical mechanisms responsible for this crossover can differ considerably 

^ ' from system to system, its ^statistical signatures follow a common pattern. This fact calls for a description in terms 

O " of Random Matrix TheoryO^B. Although this transition has already been addressed by numerous authors over the 

1^1 years, there are still many formerly unknown features being found in new areas prompting new interest and research. 

In particular, there is a constant challenge for theorists to check whether these new aspects are in agreement or in 

^-H , disagreement with the predictions of Random Matrix Theory and, whether or not they can be understood analytically. 

K*" ' Because of its statistical foundations. Random Matrix Theory is a very ambitious concept since it aims, very often 

^^ , successfully, at a unifying description of statistical features of very different physical systems. 

^^ ■ Recently, such a challenge has been posed by new developments in quantum chaology and condensed matter physics. 

*vj In studying a random matrix model describing the transition from regnlar to chaotic behavior which was motivated by a 

f^ ■ model of rotational damping in nuclear physicsQ, Persson and Aberg3 and Mizutori and AbergEI identifled, numerically 

[~^ ' and phenomenologically, a critical length scale Lc (there referred to as Lmax) in the energy spectrum. Below this scale, 

0^ , the long range correlations are chaotic and above it they become regular. This length scale was related to the spreading 

-4-^ ■ width or "energy localization length" of the wave functions. The new findings in condensed matter physics require 

Jr . a more detailed excursion. In the last ten years random matrix theory has found wide applications in the modeling 

of transport and equilibrium properties-pf mesoscopic systems. Prominent examples-include electron transport in 

quasi-one dimensional disordered wireal3 and persistent currents in mesoscopic ringglil. The general assumption in 

all these approaches was that electron-electron interactions can be neglected. In viewjofr-the considerable discrepancy 

Q I betwepi the relatively small currents calculated for rings with noninteracting electronsta'tll and the large experimental 

^ . valuealj it has become generally accepted in the last few years that interactions play a vital role in the persistent 

l^ ' current problem. Quite recently, evidence has been put forward that interactions might also change our present 

. ^H understanding of transport and localization in mesoscopic systems considerably: Shepelyansky predictedtfl that two 

^^ interacting particles in a Id disordered chain can be extended on a scale L2 far exceeding the one-particle localization 

H ] length Li. Subsequent workEj'EEl has quickly led to a definite confirmation and a more detailed understanding of this 

- - I phenomenon. As far as the two-particle effect in one dimension is concerned, there is probably only one unresolved 

issue left: the question of the parameter combination that governs the enhancement of L2 as compared to Li. 

Measuring lengths in units of the lattice spacing, and the strength U of the (Hubbard-type) interaction in units pf 

the hopping integral, Shepelyansky found L% oc U'^Lf. On the other hand, a microscopic, numerical calculationtZl 

revealed that L2 oa \U\Ll. In a, recent paperta, an attempt was made to resolve this discrepancy by studying both an 

effective Hamiltonian for the two-particle problem and a microscopic model. The effective Hamiltonian, which was 

the basis for Shepelyansky's original claim, is constructed in the following way: Diagonalize the noninteracting part of 

the two-electron problem and express the microscopic Hamiltonian in the basis of two-electron product states. The 

resulting representation consists of a diagonal contribution containing the eigenvalues of the noninteracting problem. 



and an ofFdiagonal contribution originating from the interaction operator. With the crucial assumption that both the 
diagonal and the ofFdiagonal matrix elements can be chosen to be random variables we arrive at a random matrix 
model that is precisely of the form to be studied in the present paper. The relative strength between diagonal and 
off diagonal matrix elements is determined by the interaction strength U and the one-particle localization length Li 
(i.e., the strength of the disorder). The quantity of interest in this matrix model is the energy scale Ec at which the 
spectral correlations deviate j&om pure random matrix behavior. This "spreading width" or "two-particle Thouless 
energy" was shown by ImrjO, who generalized the famous Thouless scaling picture to two-electron transport, to 
determine the two-particle localization length via the relation 

^-92{L,) = ^. (1.1) 

Here, 32 is the "two-particle conductance"Ej and A2 the two-particle level spacing. If one employs Fermi's Golden 
Rule to estimate Ec /A2 ex U^J^ (see Ref. [l5| ), Eq. (1.1) immediately leads to Shepelyansky's original relation 



L2 oc U^Lf. The main result oftEl was to point out that Fermi's Golden Rule is inappropriate as long as the off- 
diagonal matrix elements are sufficiently weak. Instead, it was found numerically that there exists a regime where 
Ec /A2 oc |C/|\/L]"- While this observation does not directly explain the numerical scaling lawE3 L2 (x \U\Li it hints 
at a possible origin of the linear dependence on U. In a recent papenij Jacquod, Shepelyansky, and Sushkov argue 
that Ec /A2 oc \U\ around the middle of the spectrum (i.e. in the vicinity of the band center) while Ec /A2 oc C/^ 
otherwise. This offers an alternative explanation for the different numerical findings. The quantity Ec /A2 should 
of course be identified with the critical length scale Lc in Rcfs. ||,||. 

It is the purpose of the present piffi^r to derive analytical formulaa_that can serve to study the critical length 
scales discussed in quantum chaologyQ'u and condensed matter physicsEj. To this end, we make use of an integral 
repijeaentation of the exact two-point correlation function for the effective Hamiltonian recently derived by one of 
uaHj'tiJ. In taking such an analytical approach we wish to present a unifying discussion which is applicable to the 
findings in both fields, in quantum chaology and condensed matter physics. Therefore, since we want to make this 
presentation self-contained and readable for a broader audience, we review and comment on earlier analytical results 
in Sec. || before we focus on the two different regimes in Sec. [11. In Sec. IV, we compute the local level density 
exactly for finite level number. We summarize and discuss our results in Sec.^. 

II. MODEL AND OUTLINE OF METHOD 



In Sec. [I A, we compare those two random matrix models describing a crossover from Poisson regular to chaotic 



fluctuations which are currently being used in the framework set in the introduction. We discuss earlier results for 



the joint probability density and the correlations in Sec. II B, In Sec. II C, we briefly review the analytical solution 
derived in Rcfs. |2^j2l|. 

A. Two Models - Same Physics 

The Hamilton matrix of our physical problem is represented by an iV x iV matrix H with random entries. To 
achieve meaningful statistical results, we have to average of the ensemble of all those random matrices which satisfy 
the physically relevant boundary conditions. Thus, the physics of such an statistical model is uniquely determined 
by choosing the probability density function of the matrices H. Here, we are interested in a crossowsr or transition 
between different regimes implying that we have to go beyond the commonly used Gaussian EnsembleaJ which describe 
fully chaotic systems. We will mainly study systems under broken time-reversal invariance since, first, some presenliy, 
investigated systems, e.g. mesoscopic rings in a magnetic fieldlij, belong to this class and, second, recent result£3'Ell 
make a detailed analytical discussion of this situation possible. For our random matrix model this means that we have 
to work with Hermitean matrices whose distribution reaches the limit of the Gaussian Unitary EnsembleEl for some 
value of the crossover or tmnsition parameter. There are two models which are currently being used in this context. 

Weinmann and Pichardll3 studied numerically the two- level correlations induced by the probability density function 



n— 1 n<m 



p(^)(if,A.) = n^exp(-iJl) n y^^exp(-(l + /.)Rc2H„„0 
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cxp (-(1 + fi)lm^Hn„i) 



(2.1) 



where we used our freedom to fix the energy scale by choosing a value of 1/2 as the variance of the Gaussian 
distributions for the diagonal elements. The parameter /i alters the variance of the distribution for the off-diagonal 
matrix elements in such a way that this ensemble interpolates between a pure Gaussian Unitary Ensemble (GUE) 
for /i = and a Poisson Ensemble for /i — > oo. In the latter case, all off-diagonal matrix elements are weighted by 
(5-functions- The parameters of this model are related to those of the physical problem of two interacting particles 
as followslla. The matrix dimension N is given by L'^{Lf -f l)/2, the number of symmetrized two-particle states 
in a d-dimensional system of size Li. The variance a^ of the diagonal elements is defined by the bandwidth B 
of the noninteracting problem, ct^ ~ B^ jZ. Normalizing this variance to 1/2 as we have done above amounts to 
measuring all energies in units of yj2/iB. In these units the variance of the off-diagonal elements is given by 
{l + ^i)-^^Q{U/BfL-^''. 

We introduce a new parameter v — 1/^/T 
form, 



/i and write the probability density function (2^) in a more convenient 
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P'^^\H,U) ^ P^^\H,u)\{p'^''\Hnn.v) 



P(l)(i7,iy) = 






(2.2) 



The distribution P^^'{H, v) is a normalized Gaussian distribution with variance v^ /2. The normalized distributions 
p^'^)[Hnn,i^) affect only the diagonal elements. The definition (2.2) includes (2.1), but is slightly more general since 
we do not specify the function p'-'*'(z, v) explicitly. The Poisson case is now recovered in the limit v 



whereas large 

values of v yield the GUE. In the sequel, we will refer to the class of models defined by Eq. (2^) as to model A. The 
/c-point correlation functions depending on k energies Xp, p = 1, . . . , fc in this model are given by 



^■'^(^i, 



^Xk.v) 



-l/d[i7]p(-^)(il,.)nimtr-l 



H 



(2.3) 



where the energies are given imaginary increments such that Xp — Xp — is. The volume element d[H] is, as usual, the 
product of all independent variables. 

The second model, which we will refer to as model B, is defined through the random Hamiltonian 



H(a) = ij(") + aH^^^ 



(2.4) 



consisting of a regular and a chaotic part, H^^> and H^^\ respectively, with an interpolating transition parameter 
a. T he matrices H^^^ belong to the GUE, i.e. the probability density function is given by P^^\H, 1) as defined in 
Eq. (2.2). The probability density function p("^(_ff(°)) is completely arbitrary. The calculation of the correlation 
functions in this model, 
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H{a) 



(2.5) 



requires an average over both matrices. However, due to the rotational invariance of the volume elements, the unitary 
matrices diagonalizing i?^^^ can be absorbed implying that, without loss of generality, H^^^^ can be chosen to be 
diagonal. Since we want to have Poisson statistics as a limiting case, we can write 



N 



p(0)(^(0)) ^ np(")(i/(0)) n ^(Rei7("„) )^(Imi/(°„)) 



(2.6) 



n=l 



where the function p'"' (z) is smooth but otherwise arbitrary. We recover the Poisson and the GUE limits for zero and 
infinite transition parameter, respectively. Model B has been studied by many authors, for a review see Refs. po|-p5[ 
Although it is intuitively obvious that Models A and B must be closely related, we feel it is worthwhile to work out 
the precise relationship between them. Gonsider model B. The form (pTq) of the correlation functions suggests to use 
the matrices H — H(a) introduced in Eq. (2.4) as new integration variables, we find 



k 

lf\xu...,Xk,a) = ^ /"d[i/(°)]p("'(i7(0)) /'d[i/]p(i)(iJ-i/(o),a) J]lmtr^^^ . (2.7) 

TI" J J , Xp — H 



Since the probability density P*^^^ is Gaussian, the diagonal and off-diagonal elements are, as far as the integration 
over H^^' is concerned, fully decoupled. In other words, the integration over the purely diagonal matrix H^^' amounts 
to nothing but a convolution of each of the functions ^'•"■'(-ffnn) with a Gaussian g{Hnn — Hnn, a) which gives a new 
probability density function 



p(°)(i/(°„)).g(i/„„ - Hi',la}dH!{'2 (2.8) 

for the matrix elements on the diagonal. Due to the fact that the total probability density function is a product 
of P^^' {H, a) and the remaining Gaussians qiK eHnm, ct) and g{ImHnm, ct) for the off-diagonal matrix elements, it is 



obviously of the form (2J) or, equivalently, (|2.2| ). Hence, we have convinced ourselves that, after a proper re-definition 
of the transition parameters, model A and model B describe the same physical situation. We may therefore confine 
our discussion to model B whose analytical properties are, advantageously, already studied in great detail. Thus, from 
now on, we drop the upper index (B). Note that these considerations can easily be generalized to the cases of real 
symmetric or quaternion self adjoint matrices. 

Finally, we would like to emphasize that there are more models which describe a transition from Poisson to GUE 
statistics. We mention two receai; works, a new formal model by Moshe et alBj and a study in the context of disordered 
systems by Altland and FuchsEZl. The fact that the two limiting cases, the Poisson and the GUE statistics, are the 
same in all these models does not imply that the interpolation between these limits is necessarily the same. The 
concrete physical situation dictates which model to choose. 

B. On the Computation of the Correlation Functions 

A quantity of central interest in every random matrix model is the joint probability density function PE{X,a) of 
the eigenvalues. In the case of Hermitean matrices, this function can be worked out explicitly for our model B. We 
diagonalize the matrix _ff by a unitary matrix U such that H = W XU where X — diag(xi, . . . , xn) is the diagonal 
matrix of the eigenvalues. To express the probability density as a function of the latter alone, we have to do the 
integrals over H'^^'^ and over [/, 

P£(A:,a)d[X] = /"d[ij(°)]pW(ij(°)) /" P''^\H-H<^"\a)d[H] . (2.9) 

J JU(N) 

The off-diagonal elements of if '^"^ are integrated out trivially leaving us with the diagonal part which we denote by X'^'^^ . 
As is well known, the volume element of H is given by d[H] = A^^{X)d[X]d^{U) where Ajv(X) = HrKml^" ~ ^™) 
is the Vandermonde determinant. Since the function P^^^ is Gaussian, the integral over the unitary group with Haar 
measure d^{U) is just the Harish-Chandra Itzykson Zuber integralEa which has been used already in Refs. HpSl to 
compute joint probability density functions of the eigenvalues. We arrive at 

PEiX,a) = -l|d[x(")]p(«)(X(°))det[ff(x„-x(^),a)]„,™.i,...,^^^^g^ (2.10) 

where g{z,a) is the normalized Gaussian with variance a^/2 which was introduced above. Note that, due to the 
symmetries of the integrand, this can be written more conveniently as 

P^(X,a) = -^|d[X(«)]p(«)(X(°))exp(-i,tr(X-X(°))^)^xW_. (2.11) 

The Gaussians yield (5-functions for vanishing variance which implies that the Poisson case is recovered in the limit 
a —f 0. For very large transition parameter a — > cxd, the GUE joint probability density function has to re-emerge. 
This limit is non-trivial but can be checked using the methods developed in Ref. |2J. 
The k-level correlation function is given as the N — k dimensional integral 

/■+00 I' + oo 

Rkixi,.. .,Xk,a) = / dxk+i--- dxN PE{X,a) (2-12) 



allowing one to use the formula (2.11) to compute them. This approach was taken by LenzlHj who managed to express 
the most interesting two- level correlations R2{xi,X2,a) as a four dimensional integral. His calculation is exact for 
every a and every N. Unfortunately, the integrand involves ratios of iV-dependent determinants which rendered the 
evaluation of the physically interesting limit of infinitely many levels impossible. Another approach was recently 
put forward by PandeycJ who chooses Dyson's Brownian motion model as a starting point. He arrives at a two 
dimensional integral representation for the two-level correlation function on the unfolded energy scale, i.e. in the limit 
N —f oo. Although his derivation is only sketchy, we have no doubts that his result is correct. Apparently, he avoids 
the explicit use of the Harish-Chandra Itzykson Zuber integral by conducting a highly non-trivial re-summation of 
certain representation functions of U {N). Yet another technique, the graded eigenvalue methodj-was used by one 



of the present authors in Refs. p0|j21| . It is a variant of Efetov's famous supersymmetry methodtj which relies on 
the observation that averages over a random potential can be most efficiently done by mapping the original field 
theory onto superspace. In Ref. pl^ it was shown that averages over random matrix ensembles lead to the same model 
in superspace. Thus, supersymmetry can be viewed as something like the "irreducible representation" of Random 
Matrix Theory. The graded eigenvalue method allows therefore for a much faster derivation of spectral correlation 
functions, even for finite level number, than any other method we are aware of. Moreover, its results go beyond 
those of Refs. P4,pq. First, higher correlations can be studied, the graded eigenvalue method yields a 2k dimensional 
integral representation for the /c-level correlation functions. Second, these integrals have a very convenient and rather 
compact form even for finite level number N. Third, the unfolding, i.e. the limit iV — > cx) is trivial and, fourth, 
the two-level correlations on the unfolded scale become a double integral for completely arbitrary initial probability 
density functions p('')(i?(")). 



C. Supersymmetry Method 



In order to make our presentation self-contained for those readers who are familiar with the main concepts of the 
supersymmetry method, we will briefly sketch the crucial steps that were taken in Refs. pQ,P^. Readers with little 
background in supersymmetry are referred to Refs. ^,^l| or, especially regarding the graded eigenvalue method, 
to Ref. p^. By including the real parts of the traces of the Green functions in Eq. (2.5), we define the functions 
Rk{xi, . . . ,Xk,a) which always allow the reconstruction of the physically interesting functions Rk{xi, . . . ,Xk,a) as 
linear combinations of the functions R^. {xi , . . . ,Xk,a). Advantageously, those can be written as the derivatives 



Rk{xi,. 



1 



Qk 



,Xk,a) 



(2-)'=np=i5J, 



Zk{x + J, a) 



(2.13) 



j=o 



of a normalized generating function Z^ {x + J, a). The energies and the source variables are ordered in the diagonal 
matrices x = diag(a;i, xi, . . . ,Xfc,a;fc) and J = diag(— Ji, -|-Ji, . . . , — J^, -|-Jfc), respectively. The physically relevant 
correlations Rk{xi, . . . ,Xk,a) are generated by the function ^Zk{x + J, a) where the symbol 3 stands for the proper 
linear combination. By using the supersymmetry method, the average over the GUE can be performed directly and 
the generating function acquires the form 



Zkix + J, a) = J d[H^°^] PJ^\H^°^)J d[a] Q,(a, a) 

detg-i ((a;± + J - a) (g) In - hk (S) i?^"^) 



(2.14) 



where cr is a 2A: x 2fc Hermitean supermatrix and In and l2k are the N x N and the 2fc x 2fc unit matrices. The 
function 



Qfe(CT,a) 



_ r,fe(fc-l) 



1 

exp I rtrgcr' 



(2.15) 



can be viewed as a normalized probability density function for the supermatrices. It is useful to shift the matrix 
X + J into this Gaussian function and then to diagonalize the supermatrix, a = u~^su by a superunitary matrix 
u where the k eigenvalues Spi, p — 1, . . . , fc in the boson boson and the k eigenvalues iSp2, P = 1, . . . , A: in the 
fermion fermioB-sector are ordered in the matrix s = diag(sii, isi2, . . . , Sfci, isk2)- The volume elements reads d[a] = 
Bl{s)d[s]d^{u)c3 where d[s] is the product of the eigenvalue differentials and d^{u) is the invariant Haar measure of 
the unitary supergroup. The square root of the Jacobian, here referred to as Bcrezinian, is given by the determinant 
Bk{s) = det [l/(spi — isq2y\„ „^i f.- In these coordinates, the generating function for the entire transition can be 
written as 



where the function 



Zk{x + J,a) = f Qk{u-hu-x-J,a)zl°\s)Bk{s)d[s]dfi{u) (2.16) 

ZJ^\x + J) = /"rf[i/(")]P^"^(i/(°))detg-i((x± + J)®l^-l2fc®i/(°)) (2.17) 

is the generating function of the correlations of the arbitr ary e nsemble defined through P'^°)(_ff("^). Obviously, the 

required average over the unitary supergroup u in Eq. ( 2.16 ) is just the supersymmetric version of the Harish- 



Chandra Itzykson Zuber integral. This observation is the main ingredient of the graded eigenvalue method. Collecting 
everything, and performing the derivatives with respect to the source variables as in Ref. ^, we find 

Rk{xi,...,Xk,a) = -^^ f Gk{s - x,a)^zl^\s) Bk{s)d[s] (2.18) 

for non-zero a. The case a = is trivial by construction. Here, the kernel is given by the Gaussian 

1 / 1 . , .2 



Gk{s~r,a) = ^^ exp ^trg(s - r)M . (2.19) 

In order to calculate the generic fluctuations, we have to unfold the correlation functions for large level number N by 
removing the dependence on the level density. We define new energies ^p = Xp/D, p = 1, . . . ,k where the mean level 
spacing D is of the order l/viV. The transition parameter a is defined on the original energy scale and has therefore 
to be unfolded, too. The new, universal transition parameter A = a/D was first introduced by PandeyO. The 
fc-level correlation functions on the unfolded scale Xk{^i, ■ ■ ■ ,^k,^) = hmAr^oo D^Rk{xi, . . . ,Xk,a) are t hen generic, 
i.e. translation invariant over the spectrum. It is useful to unfold the integration variables s in Eq. ( p.l§| ) by making 
the rescaling s -^ s/D. We arrive at 



6, A) = ^^ /Gfc(s-e,A)5zf (5)i?fe(s)4 



Xfe(ei,...,6,A) = '-^ I Gkis-^,X)^z'^>{.s)Bkis)d[.s] (2.20) 

for non-zero A where the unfolded generating function of the arbitrary correlations is given by z^, (s) = 

lirai\j^oc Z^ (Ds). Hence, we have expressed the unfolded fc-level correlation function for the transition from ar- 
bitrary to GUE fluctuations as a 2k-io\d integral. Note that the limit of infinitely many levels is trivial due to the 
application of supersymmetry. The Gaussian kernel remains unchanged when going from the original to the un- 
folded energy scale. It was shown in Ref. ^ that the generating function obeys an exact diffusion equation which is, 
remarkably, the same on both energy scales, the Gaussian kernel is the Green function of thia diffusion. 

Here, we are mainly interested in the two-level correlations. Due to translation invariancefi on the unfolded scale, 
they depend only on the difference r = ^i — ^2 of the energies. Therefore, it turns out to be useful to make the 
change si = su + S21, ti = s\\ — S21, si = S12 -I- S22 and ti = S12 — S22 of the integration variables. Since 
S5Z2 ('S) depends only on the differences of the eigenvalues by construction, it cannot depend on si -I- is2- Thus, we 
write 3z2 (s) = Szj (si — is2,ti,t2) which, for reasons of consistency, should be even in each of the differences ii 
and t2. By using a standard integral theorem of complex analysis, the integrals over si and S2 can be performed 
straightforwardly and the two-level correlations can thus be cast into the form 

-t-oc -l-oo 

^2(^,A) - ^/ y'exp(^-^(i?+t^)^sinh^sin^^^^^l^3zf(0,ti,t2)dtidi2 (2.21) 

—00 —00 

which is exact for every non-zero value of the transition parameter A. 

We emphasize that all results derived so far are correct for arbitrary initial correlations Rj, {xi,...,Xk) or 
Xj. (^1, . . . ,^/c). We now apply them to the case of the Poisson, i.e. correlation- free initial spectrum, whose probability 



density function is given in Eq. (2_^). A straightforward calculation yields 
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where R\ (x) is the Stiltjes transform of R\ (x). We point out that this is exact for all values of N. By choosing 
Rl (x) — Rl (x) to evaluate the limit N -^ oo, we find 



,(0) 



(s) = J|exp(=F«7r6p(s)) 



(2.23) 



p=i 



The signs are determined by the choice of the sign of the imaginary inc reme nt in the Green function. In the coordinates 
introduced above, the initial condition for the two- level correlations ( 2.21 ) takes the form 



Qzi°\0,h,t2) 



ti+tl 



-Re exp —ITT— 
2 V V 2^1 



- 1 



(2.24) 



which still involves an imaginary increment. In Rcfs. pO| , pT| , Eqs. (2.21) and (2.24) were used to construct a two 
dimensional integral representation for X2(r , A) in terms of Besscl functions. However, for the discussion to be 
performed in the next section, the form (2.24) is the more convenient one. 



III. TWO DIFFERENT REGIMES 



Both the numerical findingstS described in the introduction and the structure of the integral representation (2.21) 
suggest to distinguish between the two cases A ^ 1 and A ::^ 1. In these limits the exact but very cumbersome 
expression (2.21) can be substantially simplified. After introducing the proper observable in Sec. Ill A, we deal with 
the small A and large A regimes in Sees. IIIB and IIIC, respectively. 



A. Spectral long range correlations 

The spectral long range correlations are particularly well suited for the study of fluctuation properties in systems 
which undergo transitions. The level number variance S^(I/) and the spectral rigidity ^^{L) described, as functions of 
the interval length L on the unfolded scale, the fluctuations of the integrated, unfolded level density around its smooth 
value which is just L. Poisson regularity results, for both observables, in a linear behavior whereas chaotic cnriiftlations 
yield a logarithmic dependence on L. It is intuitively clear and was shown in the numerical simulationsoQclH that 
the admixture of chaotic features becomes first visible on shorter scales in the spectrum. In other words, the larger 
L, the stronger have the chaotic features to be in order to bring the long range correlations close to the logarithmic 
behavior. Thus there is a critical interval length, near which S^(i) or /S.^{L) "bend away" from the chaotic behavior 
to show a linear characteristic typical for the lack of correlations in the Poisson case. This critical length Lmax or Lc, 
respectively, was introduced and studied in Refs. g,^,^ and is also the quantity of interest in our discussion. In Ref. ||, 
the critical length imax was interpreted as the typiical spreading width of the wave functions on the unfolded energy 
scale. In the context of condensed matter physicstS, the critical length L^ was identified as the Thouless energy in 
units of the mean level spacing. Obviously, these two viewpoints describe closely related physical situations which are 
formulated in a unifying language in the random matrix model. 

We now wish to acquire information about the dependence of the pdtical length on the transition parameter A 
from the analytical discussion of the previous section. As is well knownEJ, the level number variance is related to the 
two-level correlations by the formula 



Y?{L,X)=L-2 j {L - r)Y2{r, X)dr 

Jo 
X2{r,X)^l-Y2{r,X) 



(3.1) 



where Y2{r,X) is referred to as the two-level cluster function. The "bending away" from the logarithmic behavior is 
accompanied by a change of the curvature of S^(i, A). Therefore, we may identify the point at which the curvature 
changes with the critical length. Thus, we simply have to investigate the second order derivative. 



^j:^{l,x) = -2Y2{l,x) , 



(3.2) 



which is, by construction, just the two-level cluster function. Therefore we can alternatively determine the zeros of 
Y2{L, A) or directly investigate the function E^(L, A) in order to identify the scale Lc- 



B. Small A 



The small A regime can be studied by perturbation theory. This approach has already been taken in Refs. |22|j2S 
with the result 



X2ir,X) ~ - exp 



r 
2A^/ ./o 

,i2 



r/X 



exp 1 y ) dr] 



exp 



V 



sm-dr; 



(3.3) 



which is a universal funct ion of r/X. It was checked in Ref. Ellthat the analytical expressions (2.21) and (2.24) indeed 



yield the approximation (3.3) for small A. Since Y2{r, A) is the deviation from the Poisson behavior which is just unity, 
we only have to find the zero of this function in the variable r/X. We find r/X w 1.3 implying that the critical length 
behaves as 



1.3 A 



(3.4) 



for small values of A. Thus, the critical length is indeed linear in A for small values of the transition parameter in 
agreement with the numerical investigations in Ref. |l§| . 



C. Large A 



To study the large A regime we evaluate the asymptotic ex pans ion of X^ir, A) to first order in the inverse transition 
parameter by a saddle point approximation. To bring Eqs. ( ^.2l|) and (2.24) into a form which is manifestly suitable 
for such a procedure we substitute according to ti = X^t[ and ^2 = -^^^2- Omitting the primes we obtain 



X2{r,X) = 



r3A2 



OO /' + 00 



— OO J — OO 

tlt2 



dtidt2 exp ( {ti +^2) ) sinhrti sinrt2 



C^2 ^^2\2"^^2 (0'^ *l'-^ *2) 



(3.5) 



with the initial condition given by Eq. ( ^.24 ). A direct saddle point approximation still meets with the probl em that 
there are solutions of the saddle point equations for which t\^ t\ — Q. To avoid the ensuing singularity in ( p.q ) we 
introduce an additional parameter /3 and a function 



^2(r,A,/3) = 



r3A2 



00 /'+00 



—00 -' —00 



dtidt2 exp ( f3{il + ^3) ) sinhrti sinrt2 



,(0) 



ht2^Z^"'{0,X'h,XH2) , 



(3.6) 



in terms of which we have 



X2(r,A) = — lim / d/3' / d/3X2(r, A,/3) 
4 v^^Ji J pi 



(3.7) 



In the integrand of Eq. (3.6), the product sinhrii sinrt2 can be replaced by the expression — iexp(— rii) exp(— iri2) 



because the additional contributions vanish upon integration (they make the integrand odd in either ti or ^2 or both) . 
We obtain 



^2(r,A,/3) 



7r3A2 



CO /' + 00 



00 J —00 



dtidt2 exp ( — TrP{ti + ^2) ^ ''^i ~ ^'"'■2 



t + 1 
tit2 ( exp ( -i-kX^— — 3^ 



- 1 + c.c. 



(3.8) 



The contributions originating from the exponential and the constant in the square brackets, respectively, have to be 
treated separately. We begin with the more complicated term, namely 



xt\rA,P) 



7r3A2 
A2 



+ OC p+oc 



dtidt2 exp(-£(ti,t2)) tit2 

j-2 I 7-2 



2*1+^2 



>C(ti,t2) = Y/3(i?+ti)+i7rA ^^_ 



rii 



irt2 



(3.9) 



To determine the saddle points we can neglect the imaginary increment for ti in the second term of C(ti,t2). The 
saddle point equations dt^C = and dt^C ^ Q read 



X^pti 






>2^2 



+ r = 



y(it2+iTTV—+ir = Q . 



With the abbreviations 

7 = iTrA 
these equations can be cast into the form 



and 



t = X^(3t2 + ir 



(3.10) 



(3.11) 



7^2 

t 

3 I /_,2 o„_A+ o„-„-,2 



^ + (7 - 2r'7)t - 2ir7^ = 
The roots of the cubic equation are given by 



(3.12) 



tW- 



t(2) = z|(i + v/T^^^) , 



(3.13) 



.(«) 



in terms of which the solutions tl 2 (* = 1, 2, 3) for the original variables are easily recovered, 
jroximation proceeds according to the following sche 

X^^\r, A, /?) « -^5. ^ t^^hi^^ exp (-£(4^'\ 4^'^)) 



VdetAU) 



(3.14) 



where 



A(^") = 



d^C uU) i-U)\ d^C(Aj) j-{j)\ 

dt2dti\''l '''2 I 9tfv''l ''2 J 



(3.15) 



is the matrix defining the form of the Gaussian fluctuations around the saddle point. The ap prox imation (3.14) is 
valid for A S> 1 and all r . Expressing the exponent, the pre-exponential terms, and detA in ( 3.14 ) in terms of the 
saddle point values (3.13) we get for j — 1,2,3 



/:(4'\4''^) = 



2A2/3(iO)): 



.[^-itU)f{r + iM) , 



A3)Aj) 

'•1 ''2 



yiO) 



A4/32 
detA(^' = {X^pf (l 



i^-'w) 



tU) 



iO") 



(t(J))2 



(3.16) 



At this point we notice that i^^-* is not a valid saddle point. This can be most easily seen in the extreme GUE limit 
(A — > cx), I7I — * 00), where i'^^-' — > 2ir. Obviously, detA^'^) < in this case, i.e. the matrix A^^^^ is not positive definite. 
Therefore, t^^' does not correspond to a maximum of the integrand and will be discarded henceforth. 
Before we proceed we have to investigate the second contribution in (p^). 



XT{r.\fi) 



C.{hM) 



7r3A2 



+ 00 z' + OO 



dtidt2 exp 



00 -/ — CXD 



(-Aii,t2)) 



iit 



l'-2 



A^ 



/3(tf + ti) + rti+irt2 



Here, the saddle point equations are trivial and have the single solution 

r ^ ^ ir 



U =-- 



A2/3 



and 



t2 = 



A2/3 



(3.17) 



(3.18) 



Therefore, we get in analogy to ( 3.16 ) 



C{hM)^0 , 



hu^ 



(A2/3)2 ' 
detA = [X^pf . 



Finally, collecting our results from ( 3.13 ), ( p.lq ), and ( 3.19 ) and introducing a new parameter through 

T . r 

— = — ^ — = —IK , 

7 ttA^ 

we arrive after some algebra at 

l(r,A,/3)=X('^)(r,A,/?)+lW(r,A,/3) 

4 / (2iK+l + p)(3 + p) /iTTr{2iK+l + p){3 + p)^\ 

1 ^^ cxp I -TTT r:— — T^^ I + c.c. 



(3.19) 



(3.20) 



A4/33 



8VP 



4/3 



(1 + pF 



with 



p = \J\ + 8lK . 



(3.21) 



(3.22) 



The integration over the auxiliary parameter /3 (see (3.7)) is straightforward and leads to the final result of this 
section, 



X2(r,A) = l + 



(1 + Pr 



(7rr)2 (2m+l + p)(3 + /9)3^ 

(2iK+H-p)(3 + p)2 



exp iTxr 



4(l + p)2 



. (2iK+l + p)(3 + p)2 _, 

-iTxr -7; 1 

4(l + p)2 



(3.23) 



Although the decisive parameter governing the deviation from the GUE limit is k = r jiO^ , closer inspection reveals 
that the critical length L^. is still proportional to A. We thank Y. Fyodorov for this insightCJ. In a first order expansion 
in K <C 1, one finds 



X2(r,A)«l + 



1 



-2'7rr 



2(7rr)2 



^ cos 27rr + 27rrK — 1 ) 



For K ^ this reduces to the GUE limit 

X2(r,A^oo) = l 



(3.24) 



(3.25) 



For finite but small k, Eq. (3.24) shows that the oscillations in Xg.(r, A) are exponentially damped on the scale 
nr oc r^ j}? so that Lc must behave like A. Numerical evaluationsEj of formula (3.23), which did not rely on the 
assumption k <C 1, confirmed this view. At this point, an important caveat is in order. In Random Matrix Theory, 
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the unfolding procedure usually involves the crucial assumption that the transition parameter a scales like the mean 
level spacing as the level number N increases. This is formally necessary to keep the ratio A = a/Z?, i.e. the transition 
parameter on the unfolded energy scale, fixed when the limit A^ — > oo is taken. Physically, this is a well justified 
assumption which means that the chaotic admixture, although affecting the fluctuation properties considerably, does 
not change the mean level spacing. The results of Refs. E0,G1^ which are the starting point of our present considerations 
were derived under this assumption. It is not completely straightforward to relate our parameter A to the transition 
parameters of the recent studies in Refs. [l8Hl9| . Further research along t hese l ines is under wayE3. 

It is interesting to notea that under the assumption r <C A ^ A^ Eq. (3.24) can be further simplified to give 



^2(r-,A>r) = 1- 



smTrr 



7r2A2 



■ sin 7rr 



(3.26) 



Remarkably, this expression coincides exactly with a result derived by Kravtsov and MirlintZI for the spectral two-point 
function of a quasi-one dimensional disordered wire close to the universal Wigner-Dyson limit, provided A is identified 
(via Thouless' argument) with the conductance. The third term in Eq. (3.26) represents the first nonperturbative 
correction to Wigner-Dyson statistics. It is intriguing to further explore this analogy. 

Finally, we would like to mestion that we have tested the quantitative applicability of Eq. ( 3.23|) by means of 
extensive numerical simulationg23. We found that the number variance Y?[V) derived from formula (3.23) describes 
the data up to the largest spectral ranges of about L rj 140 we investigated, as long as A 3> 1. We have therefore 
confidence in the usefulness and correctness of Eq. 



(3.23) 



IV. LOCAL LEVEL DENSITY 



The local density describes the influence of an admixture on the ?Ti-th level, say, of a given spectrum. Consider the 
Hamiltonian H{a) — iJ'^"-' + aH^-^^^ defined in Eq. (2.4) with a diagonal regular part H^^\ If we avoid averaging over 
iJ^°'', we can study how the density around the rn-th level of if'"' is affected as a (and hence the chaotic features) 
increase. This quantity was studied numerically in Refs. Esl and analytically in Refs. p^,p9|. However, this analytical 
derivation relies on the saddlepoint approximation and is therefore only valid for rather large values of a, i.e. when the 
level density is sufficiently smooth. It is the main purpose of this section to show that a direct and exact evaluation, 
valid for all values of a, of the local level density is easily possible by the methods outlined in Sec. |l[ even for finite 
level numbers. With this result we can study the local level density in the regime of small values a which is not 
accessible to the saddle point approximation. Apart from this, the calculation which follows is also of general interest 
for nratrix models. In this context, we mention related work on Wigner random band matricealj, which also focuses 
on the local density of states. 

We express the local level density including its real part as the derivative 



Rim{x,a) 



1 9Zi(x, J, a) 



2tt 



dJ„ 



j=o 



of the generating function 

Zi{x,J,a) = / d[c7]Qi{a,a)detg~^ ({x^ - ct) «) Ijv - I2 «> -ff 



(0) 



J 



(4.1) 



(4.2) 



depending on a source field J = diag(Ji, . . . , Jn) which resolves the individual levels n = 1, . . . ,N. We introduced the 
matrix t — diag(— 1, +1). As usual, we diagonalize the Hernfitean supermatrix, a = u^^su, where s = diag(si,is2)- 
For the unitary supermatrix, we use the explicit parameterization 



l + /3/3*/2 /3 

13* 1 + (3* (3/2 



(4.3) 



in terms of the complex anticommuting angle f3 which allows us to cast the superdeterminant in Eq. ( |4.2| ) to first 
order in J into the form 



nX ^ 'A2 ~ ^nn ■^^ J 



„=i X- Si-H, 



(0) 

nn n=l 



1 + 2f3P* 
X ^~ Si — ilnri 



1 + 2P*P 



(4.4) 
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The derivative with respect to Jm can now be performed easily. To integrate over the eigenvalues and angles, we have 
to take the boundary contributions, pften called Efetov Wegner terms, into account. As shown in Ref. ^, this can be 
done by applying Rothstein's theoryE^I which yields the volume element 

^ ^El^dpdp* (l - P*P{s, - ^s,) f A _ ,1\\ (4.5) 

[S1-1S2Y V \osi 0S2JJ 

in which the term containing the derivatives takes care of the boundary contributions. Collecting everything, we 
arrive after a straightforward calculation at the expression 

i?i™(x, a) = - _ _ (0) + TTTir ^^'^' ") (^-^^ 

for non-zero a. The first part, the Efetov Wegner term, is simply the unperturbed Green function of the m-th level. 
The second part can be written as the derivative of the function 



Using the identity 



N N 

-Q X - vn = 1 + ^ ''" ~ ^" TT ''' ~ ^" (4.8) 

^^^ a; - ^„ ^^^ /^„ - a; ^_^^^ ^ - Ai« 



for ^„ = si + Hnn and i^n = is2 + i/nn j we obtain 



^ "+- '■+°o , , / . X 6„(s,iJ(")) 



6„(.,ff(o))=nfi 



S 






H^ + T 



rr(0) _ rr(0) 



(4.9) 



This can be worked out further by introducing the permutation invariant symmetric functions ci{h), I = 0, . . . , N — 1 
of a set of variables hi, I = 1, . . . ,N — 1 which are defined through 



JV-l N-l 



II (1 + hiz) = J2 ci{h)z' . (4.10) 



1=1 1=0 

In particular we have 

TV-l AT-l 



Co{h) = l, ci(/i)=^/i/, C2(/i) =^/i;/i;', ••■, CN-i{h) ^ Yl hi . (4.11) 



Kl' 1=1 

An) 



\/{H\i — Hnn) for I ^ n, we can express the first of Eqs. ( |4.9| ) in the form 



Since 6„(s, H^^'') has the structure of the generating function on the left hand side of Eq. ( 4.10| ) if we choose h^ — 
for I ^ 71, we can express the first 

N N-l 



n=l 1=0 

-00 /' + 00 / -] 

„2 I „2 



^ / 2 , 2^ \ 1*1 ^ **2J 

(0) 



— 00 J —00 



dsids2exp( -{si + si)]^ ^-^ . (4.12) 

X - Si- Hn 



^nn 



The integral over S2 yields precisely the Hermitean polynomials 



Hi^i{z)^^ / exp(-C')(^-<)'"'dC (4.13) 
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where, in particular, we have -ff_i(z) = — \/7rexp(z^)(l + erf(z)). The remaining si integral can be viewed as the 
Stiltjes transform 



A;_i(z,a) = - / dsi 

TT J-oo Z - Si 



Hence, both integrals can be done and we arrive at 

N N-l , ,w , 



V "" „=i /=o 



(4.14) 



(4.15) 



We are mainly interested in the local level density, i.e. in the imaginary part Rim{x,a) — Imi?im(x, a), which is 
according to Eq. (4.6) given by 



Ri^{x,a)^S{x~H^^l) + 



d 



uiimrr 



■W{x,a) 



(4.16) 



where W{x, a) = linW{x, a). Since the imaginary part of the Stiltjes transform is here just an integration over a S 
function, we find 



N N-l 



V „=1 1=0 



^^ Q(/^(")) exp f- ^""^f"^' ) H,^i ((. - H^^2)/c 



(4.17) 



for an arbitrary matrix H^^K There is a subtle point here to be remarked. According to Eq. (4.12), the function 



W{x, a) is ev en in a. This, however, is not immediately obvious from Eq. (4.17). When we introduced the Hermitean 



polynomials (4.13), we used 82/0 as the new integration variable which eventually led to the term a in Eq. (4.17). To 
be mathematically cleaner, we should have used S2/|a| which would have produced the term |a|'. Hence, we always 
assume a > 0. 



The form (4.17) of the function W{x, a) can be viewed as an expansion in the increasing complexity of the contri- 
butions due to the matrix i/'"^ which is reflected in the term a''ci{h^'^'>). To illustrate this, we c onsider the case of 
very small a in which only the I — term contributes significantly. Hence we find with Eq. ( 1.11 ) 



1 ^ / 

W{x, a) ~ — E exp ( - 

^ n—1 \ 



^"~J^"^V -^(("-^^»^)A 



N 

.E(l + crf((x-i7W)/a) 



ra=l 



which immediately implies 



i?i„(x, a) 



exp 



(x-H, 



(0) \2' 

imn ) 



(4.18) 



(4.19) 



Thus, for a very small admixture of H'^^\ the states of H'^^'^ are smeared out with a Gaussian shape. T his ca n only 
be true as long as these Gaussians do not overlap. As the admixture is increased with a, the Gaussian (4.19) of the 
?TT--th level slowly starts feeling the influence of the other levels of i/'-^-' and higher and higher orders in I have to be 
taken into account. The complexity of this interaction is described through the more and more entangled structure of 
the symmetric functions ci{h^^^) as I increases. We mention in passing that, as is easily shown, the odd contributions 
in / can be neglected, provided, the m-th level lies in the middle of a very lon g s p ectrum. The shape of the local level 
density changes due to this coupling to the other levels of iJ^"-* . In Ref. ^,|l8|, it was shown numerica lly th at this 
shape is Lorentzian, even for relatively small admixture. This is no contradiction to the Gaussian shape ( 4.19| ) which 
only applies in the case of very small a. Pictorially speaking, the coupling to the other levels of i/*^°-' lifts the tails 
of this Gaussian which then becomes a Lorentzian. In Ref. ^ this Lorentzian shape was analytically derived in the 
case of strong admixture. 

A further insight into the expansions ( 4.15| ) and ( [4.17|) is provided by the following consideration. Take a harmonic 
oscillator with mean level spacing D^^> as initial condition H^^' , i.e. a picket fence spectrum with equidistant levels. 
The symmetric function ci{h^^') scales now as 1/{D^^')^ implying that the expansion parameter is now given by 
Q,^£)(0) ^ This is a measure of the chaos producing interaction on the scale of the mean level spacing of the unperturbed 
spectrum. Qualitatively, this does also apply to arbitrary initial conditions H^^\ 
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V. SUMMARY AND DISCUSSION 

In this paper we have investigated spectral properties of a particular class of random Hamiltonians, namely the 
superposition of a diagonal matrix with Poisson statistics and a GUE. The motivation for this investigation came, 
apart from the obvious relevance for quantum chaology, from the fact that this class of matrices can serve as an 
ensemble of effective Hamiltonians for problems of current interest in both nuclear and condensed matter physics. 
After showing the equivalence of two alternative definitions of this random matrix ensemble, we analytically discussed 
two regimes, small values of the parametei; A, corresponding to small values of the electron-electron interaction U in 
the problem of coherent pair propagatiorO, and the regime of large values of A. In the limit A ^ 1 wc derived a novel 
explicit expression for the two point correlation function. In the crossover regime A « 1 this expression and the already 
known perturbative result for small A differ from each other and do no longer give a quantitatively correct picture. 



Comparison with the exact (albeit difficult to handle) expression (2.21) however reveals that the qualitative features of 
the crossover are still well described. This opens the possibility to study analytically and in detail statistical measures 
like the miniber variance S^. Such a discussion, furnished with extensive numerical simulations, will be presented 
elsewhereEj. 

In Ref. @, the transition from regular to chaotic fluctuations was studied using a random matrix model in which 
the total level density undergoes a transition from a sharp Gaussian, in the regular case, to the Wigner semicircle. 
With the help of Fermi's Golden Rule, the critical length L^ax was estimated as a function of the model parameters. 
The result in Ref || cannot be directly related to our discussion of the critical scale Lc- The chaoticity parameter A in 
Ref. g is defined on the original (and not the unfolded) energy scale, calling for a proper rescalings of the numerical 
and the analytical results. In Ref. an additional time dependence was introduced to study energy dissipation. 
Depending on the value of the chaoticity parameter A a time scale i* is identified. It distinguishes between a regime 
of anomalous diffusion proportional to t^ for t < i* and normal diffusion for f > t^; of the energy in the eigenbasis at 
t — 0. The precise relation of these interesting findings to the results presented here has yet to be investigated. It is 
quite plausible, however, that the different (chaotic and regular) regimes in the spectral statistics manifest themselves 
in this difference in the energy diffusion. 

Finally we have shown how to calculate, for finite TV, arbitrary values of a, and without approximation, the local 
density of states of the matrix ensemble. For extremely small values of a the shape of the local density of states 
turned out to be Gaussian. For larger Oj^ yakift s, when the levels on the diagonal of the matrix are no longer strictly 
isolated, a Lorentzian shape is expect edtaHH. 

We believe to have demonstrated in this paper that many statistical properties of the ensemble of pajadom matrices 
studied here are well under control. This is particularly important in view of recent devclopmentsliSula, where the 
ensemble considered in our work plays a central role. 
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